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ABSTRACT 

We present a new analysis of the dust obscuration in starburst galaxies at low and high redshift. This study 
is motivated by our unique sample of the most extreme UV-selected starburst galaxies in the nearby universe 
(z < 0.3), found to be good analogs of high-redshift Lyman Break Galaxies (LBGs) in most of their physical 
properties. We find that the dust properties of the Lyman Break Analogs (LBAs) are consistent with the relation 
derived previously by Meurer et al. (M99) that is commonly used to dust-correct star formation rate measure- 
ments at a very wide range of redshifts. We directly compare our results with high redshift samples (LBGs, 
"BzK", and sub-mm galaxies at z ~ 2 — 3) having IR data either from Spitzer or Herschel. The attenuation 
in typical LBGs at z ~ 2 — 3 and LBAs is very similar. Because LBAs are much better analogs to LBGs 
compared to previous local star-forming samples, including M99, the practice of dust-correcting the SFRs of 
high redshift galaxies based on the local calibration is now placed on a much more solid ground. We illustrate 
the importance of this result by showing how the locally calibrated relation between UV measurements and 
extinction is used to estimate the integrated, dust-corrected star formation rate density at z ~ 2 — 6. 
Subject headings: galaxies: starburst — galaxies: peculiar — dust, extinction — galaxies: high-redshift 



1. INTRODUCTION 

The ultraviolet (UV) and far-infrared (far-IR) emission of 
star-forming galaxies (SFGs) offers one of the most direct 
estimates of their star formation rate (SFR). The fraction of 
light emerging in the IR relative to that in the UV is mod- 
ulated primarily by t he amount of dust seen by young stars 
dMeurer et all 119991 ICalzetfi et alJ 120001 iBuat et all 120021: 
iBurgarella et al.ll2005USeibert et alJ2005HTrever et alJ2007i) . 
with seco ndary effects due to, e.g., star formation history and 
geometry (ICharlot & FaifeOOOtlKong et al.ll200ltllnoue et al] 



2006; Boissier et al. 2007; Johnson et al. 2007; Salim et al. 
2007[ IPanuzzo et alJE007[ ICortese et alj|2008t iBoselli et all 
2010). For a relatively young galaxy forming stars at a 
stable rate, the SFR is prop ortional to the UV luminosity 
dLeitherer & Heckmanl [19 95). while the degree of redden- 
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ing of the UV continuum is related to the amount of dust. 
This reddening can be expressed in terms of the UV slope, (3 
with f\ oc A' 3 , an d the total attenuati on is given by the ratio 
IRX=L IR /Luv dMeureret all 19991 M99). The "IRX-/3" 
relation of M99 is based on a large sample of local starburst 
galaxies, and has proven to be an extremely useful tool for 
estimating total SFRs from UV-only data. For the relation 
to work, it is implied that the dust must be near the UV 
sou rces and have some kind of shell- or screen-like geome- 
try dGordon et alJl2000h . 

The IRX relation is particularly important at high red- 
shift. SFR estimates based on IR, X-ray, or rad io data are 
typically available only for th e brightest objects dPope et al] 
120061 ISiana eTaT]l2008l [20091 iReddv et aljl2 006. 201 0.R10I). 
or statistically through stacks (ISeibert et al|2002HCarilli et alJ 

~et all 



2008: lReddv eTalEoTa RlO: Magdis et al. 2010a b; Ho c c. 
bold iKurczvnski et all 120 10; Rigopoul ou et alj|2010h . Sur- 



veys in the rest-frame UV can simultaneously give Lpjjy 
and j3 and hence the dust-corrected SFR ba sed on the local 
IRX relation. Using this technique, iBouwens et alJ (120091 
B09) estimated the z ~ 3 — 6 SFR density from large sam- 
ples of LBGs. They found that most of the energy output at 
z = 2.5 - 4 (z = 4-6) occurs in the IR (UV), and that 
the SFR density can be largely recovered by dust-correcting 
the UV measurements thus demonstrating the importance of 
UV surveys. However, because local calibrations of the IRX 
relation are based on galaxies that are very different from 
typical UV-selected starburst galaxies at z > 3, it is not 
clear whether they are valid at high redshift. In this Let- 
ter, we re-investigate the IRX-/3 relation for the first time 
using a unique population of nearby (z < 0.3) starbursts 
that has been shown to be similar to LBGs in most of their 
basic physical properties. These "Lyman Break Analogs" 
(LBAs) are similar in mass, age, size, metallicity, optical ex- 
tinction, and SF R (Heckman et al. 20051: iHoopes et alJl2007l 
iBasu-Zvch et alj | 2007l). have similar com pact and clumpy 
morphologies dOverzier et all 120081 120101) . similar kinemat- 
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ics dBasu-Zvch et al.ll2009l : lGoncalves et alj2010ft. and a sim- 
ilar fe edback-dominated interstellar medium dOverzier et alj 
120091 O09). The structure of this Letter is as follows. We first 
present our data and measurements (§2). We compare LBAs 
with low and high redshift SFGs having good IR data from 
Spitzer or Herschel (§3), and we discuss the possible implica- 
tions of our results on the SFR density at high redshift (§4). 

2. DATA AND METHODS 
2.1. Lyman Break Analogs UV and IR data 

The sample of 3 1 LBAs from O09 forms the basis of this 
paper. FUV luminosities, Lpuv = ALa with A = 1600A, 
were calculated from our F150LP (A C «1614A) data taken 
with the Hubble Space Telescope (H ST) Advanced Camera 
for Surveys (see lOverzier et al.ll2008l O09, for details). The 
slope of the UV continuum, Pgalex, was calculated^ from 
Galaxy Evolution Explorer (GALEX) General Release 6 data 
in the FUV (A C «1530A) and NUV (A C «2315A) having total 
exposure times of ~100-13,000 s. All LBAs are unresolved 
in the GALEX images, and we measured total magnitudes in 
a 9" radius aperture. We perform a small if-correction to ob- 
tain Lpuv an d fie alex at z = using a set of starburst 
templates. In this paper we work exclusively with the photo- 
metric measure of f3 (as opposed to the "true" f3 obtained from 
a power-law fit to a UV spectrum), that is most appropriate 
for comparing with samples for which spectra are usually not 
available in large numbers. 

We use IR data obtained with the Infrared Array Cam- 
era (IRAC), the Multi-band Imaging Photometer (MIPS), and 
the InfraRed Spectrograph (IRS) on the Spitzer Space Tele- 
scope (O09, Armus et al. in prep.). At 24/im and shorter 
wavelengths, point source photometry was performed on the 
post-Basic Calibration Data ("pBCD"). At 70 and 160/im, 
flux densities were obtained by aperture photometry on the 
filtered pBCD images. Upper limits were determined from 
the standard deviation images. The IRS spectra covering ob- 
served wavelengths in the range 5-35/j,m were normalized 
to the MIPS 24pm data. We estimate IR luminosities, Lm 
(3-1000Atm), by fitting our da ta with the model library of 
ISiebenmorgen & Kriigell (120071) (SK07). This method allows 
us to fit models covering a wide range of physical parameters, 
thus obtaining a good sense of the range of models allowed 
by our data within the measurement errors. We simulate our 
IRAC+MIPS+IRS data set (including measurement errors and 
upper limits) and calculate Lj r and its error by taking the me- 
dian and standard deviation of all the templates that best fit the 
data in the monte carlo simulatiorQ Four objects lacking a 
sufficient number of IR data points were removed from the 
sample, as well as one object with an (obscured) active nu- 
cleus identified in recent optical spectroscopy. The UV and 
IR measurements for the 26 remaining objects are given in 
Table Q] 

2.2. Redetermination of the Meurer et al. (1999) relation 
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14 While our results obtained using SK07 are consistent with those ob- 
tained from black-body fitting, we note that the SK07 library generally per- 
forms bett er for LBAs compared to the more commonly used empirical li- 
brary from Chary & Elbaz 12001). The latter often failed to simultaneously 
fit the mid-IR and the far-IR dust emission, presumably due to the fact that 
the LBAs show a greater spectral variation than typical iR-selected starbursts 
in the local universe. 



In order to be able to compare the results from differ- 
ent samples at low and high redshift in a consistent manner, 
we have re-measured the original M99 IRX-/3 relation using 
GALEX data and an updated estimate of L\r. 47 galaxies of 
the M99 sample are covered by GALEX. Details on the pho- 
tometry will be provide d in a forthcoming paper. In brief, we 
follow the procedures in Wang et ail (120101) carefully masking 
neighbours and image artefacts, perform background subtrac- 
tion, and photometry on PSF-matched images using elliptical 
apertures out to a maximum radius, r max , of 2.5 — 5xrfc ro „. 
Because r max is in most cases much larger than the Inter- 
national UV Explorer (IUE) aperture of 10" x 20" used by 
M99, we also measure the inner UV colors and fluxes within a 
maximum radius of 10" for comparison with earlier work. We 
compute Lm using the Infared Astronomical Satellite (IRA S) 
fluxes at 12, 25, 60, and lOO^rn dSanders & Mirabelll 19961) . 

3. RESULTS 

3.1. The IRX-/3 Relation of LBAs versus M99 

In Fig. Q]we show the IRX-/3 diagram for LBAs. Panels 
show the corresponding HST morphologies from O09, and 
the distribution of U V sl opes (top panel). We also show the 
values measured in i]2.2| for the M99 galaxies using an aper- 
ture radius of 10" most similar to M99, and within R ma x en- 
closing the entire source. We follow M99 and write IRX in 
terms of Apuv, the attenuation at 1600A, and the bolometric 
corrections to the total light emitted by stars (BCfuv,*) and 
dust (BCdust)- 



log 10 (IRX) = log 10 (10 



0AA FUV 



l) + lo gl 



BC 



FUV,* 



BC 
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(1) 

with BC FUV ,* ~ 1-68 (M99, Seibertet all 120051) . and 
BCdust ~ 1 for our estimate of Lm. We then derive 
the best-fit relation by performing a linear fit of the form 
Afuv = Co + Gi/3 to the data in Fig. [T] Similar to 
M99 we exclude galaxies having radii of >2'. We find 
A FUV = 4.54 + 2.07/3±0.4(IRX M 99 inner, dotted line), and 
A FUV = 3.85 + 1.96/3±0.4(IRXM99,'totai, dashed line). Our 
redetermination of IRXM99, inner is very similar to the origi- 
nal M99 relation (IRXm99.o> thin solid line). However, using 
the new larger aperture we find ~ 2 x more flux in the FUV 
for nearly all sources, and hence a lower IRX. The smaller 
IUE aperture used by M99 may thus have missed about half 
the light. On the other hand, in both cases we use a sin- 
gle integrated IR luminosity of unknown spatial distribution. 
For a subset of 12 galaxies from M99, we were able to com- 
pare the IRAS measurements with more recent determinations 
of Ljb bas ed on Spitzer photome try encompassing the en- 
tire source (Engelbr acht et al.ll2008l) . The results are in good 
agreement. If a significant fraction of the IR is due to heating 
by stars further out, using IRXm99,o rather than IRX^/gg, total 
would lead to an overestimate of the total attenuation. Alter- 
natively, IRXM99, total could be affected by the (redder) UV 
light from somewhat older stellar populations on the outskirts 
of the galaxies. In this case, the interpretation of IRXA/gg jtotQ j 
would be less straightforw ard as it requires knowledge of the 
star formation history (see lKong et a l. 2004). 

How do LBAs compare to these IRX relations? A fit to 
LBAs gives A FUV = 4.01 + 1.81/3 (implied FUV attenua- 
tions of 0-3 mag). The dispersion found in IKXlba is ~0.6. 
As shown in Fig. [T] the LBAs are thus in much better agree- 
ment, on average, with M99 (inner) than with M99 (total), but 
we note that this only applies to the range of— 2.5 < (3 < —1 



Dust Attenuation in UV-Selected Starbursts at High Redshift and their Local Counterparts 



3 



-2.5 -2.0 -1.5 -1.0 -0.5 0.0 



LB As 



1 

o 

IT 

>< 

o 






% 10 


.20 












ye 




| . 7 






% 6 






5 


. 15 




4 4 


0* 




3 






% 2 


# 12 






i 11 



-2.5 -2.0 -1.5 -1.0 

Pgalex,z=o 



-0.5 0.0 



161 ^ 26 



Fig. 1 . — : Main panel: The relation between the UV slope and "IRX" for LBAs (large symbols). LBAs referred to as "Dominant 
Compact Objects" (DCOs, O09) are marked in red, and HST morphologies from O09 are shown on the right. Open and filled 
squares indicate M99 galaxies re-measured using GALEX data (§2.2) using a small aperture similar to the original IUE measure- 
ments (open squares), and using an aperture enclosing the entire galaxy (filled squares). Fits to the data are indicated (thick black 
line: IRXlba', grey dotted line: IRX a/99, inner', grey dashed line: IRX m 99, total)- The original relation from M99, IRXm99,o, is 
shown for comparison (thin black line). Top panel: Distribution in /3 z= q for LBAs. 



probed by these samples. We conclude that the dust proper- 
ties of LBAs are most similar to those found only in the inner 
starburst cores of the M99 galaxies. 

3.2. The Attenuation of Starbursts at Low and High Redshift 

With good estimates of Lib, and Lpuv and thus the bolo- 
metric luminosity (Li, i=Lfuv + Lib), we can compare 
LBAs to other local and high redshift SFGs. This is shown 
in Fig. |2 In the left panel we show L^i versus the at- 
tenuation (Lib,/ Lfuv) f° r LBAs, ordinary SFGs, M99 star- 
bursts, and Luminous IR Galaxies (LIRGs) all at low red- 
shift, and "BM/BX" LBGs at z ~ 2 from R10. Typical 
SFGs in the local universe lie along a broad sequence in which 
the attenuation roughly follows L DO i (indicated by the dotted 
line). LBAs and LBGs follow a similar relation, but one that 
is offset by 1-2 dex toward lower attenuations at the same 
Lboi (solid line shows the best-fit relation from R10). Con- 
versely, at fixed L}, i LBAs and LBGs reach much higher 
FUV luminosities than typical local star-forming galaxies (i.e. 



'FUV 



> 0.3L* =3 , right panel of Fig. |2). 



Because LBAs appear so similar compared to z ~ 2 LBGs 
in terms of their Lboi and attenuation (Fig. |2}, it is interest- 
ing to see whether the IRX relation at low and high redshift 
is similar as well. In Fig. [3] we show a compilation of low 



and high redshift data for which accurate measurements are 
available. We first compare LBAs with the sample of z < 0.2 
SFGs detected with GALEX and Herschel from Buat et al. 
(2010), and with LIRGs from iHowell et al.l d2010l) . shown 
left. Aperture effects either in the UV or the IR are most 
likely much less of a problem compared to the M99 sample, 
due to the higher average redshifts of the two samples com- 
pared to that of M99. Therefore we will use IRXj\/9g totQ / 
as our reference (dashed line). LIRGs tend to lie above this 
relation, while the galaxies from iBuat et al.l (120101) lie below 
it (on average). For these samples, we would therefore tend 
to, respectively, under- and overestimate the attenuation when 
using IRX M 99. total, while IRX^.mner performs somew hat 
better for LIRGs (see lBuat et al.ll2010tlHowell et alj|2010l for 
detailed discussion). 

In the right panel of Fig. [3]we present an overview of the sit- 
uation at z > 2, showing LBGs, BzK a nd sub-mm galaxies 
(SMGs) at z = 2 - 3 dReddv et alj|200fl R10). le nsed LBGs 
("cB58" and the cosmic "Eve":lSiana et al][2008l [2009b. and 
stacked LBGs at z ~ 3 ( Magdis et alj|2010allbl) . LBAs and 
z ~ 2 LBGs occupy a very similar region in this IRX-/3 di- 
agram, confirming the similarity between the two samples. 
We should note that the IRX estimate of R10 is not based 
on a direct measurement of the IR emission (it is based on 




Fig. 2. — : Li, i versus the attenuation (Lju/Lpuv, left panel) and Lpuv (right panel) for LBAs (l arge circles), z ~ 2 LBGs 
(blue triangles, from RIO), M99 ( grey square s , with R < R ma x), LIRGs (open magenta circles, How ell et alj (120101) ). and 
ordinary SFGs (small circles, see Hua ng et"aT] (120091) ). The local samples (M99, SFGs and LIRGs) lie along a broad relation 
between Lf, c ; and attenuation (the result of a simple linear fit to guide the eye is indicated by the dotted line). The most extreme 
UV-selected samples at low (LBAs) and high redshift (LBGs) also lie along such a relation, but one that is offset by 1-2 dex 
towards lower attenuations with respect to the former. The mean relation at z ~ 2 from RIO is indicated (solid: BM/BX sample; 
dashed: extrapolation to intrinsically faint objects). The fraction of L{, ; emitted in the FUV is much higher for LBGs and LBAs 
compared to typical local SFGs (right panel). 



a combination of Lpuv, Lna and Lg), but statistical detec- 
tions in the X-ra ys appear consisten t with the extrapolated es- 
timates of Lib, (Sei bert et al.l l2002. RIO). We can also com- 
pare our results with the statistical detection of z ~ 3 LBGs 
in stacks at 100/im and 160/im from recent Herschel observa- 
tions. This result is indicated by the upper magenta triangle, 
which falls right in the middle of the distribution of LBAs 
and z ~ 2 LBGs. This suggests that even for the most IR- 
luminous LBGs (median Ljr of 1.6 x 1O 12 L , i.e. ULIRGs) 
the attenuation is exactly as expected based on the locally de- 
termined IRX (inner) relation. 

Direct measurements of the mid- or far-IR emission from 
LBGs are available in only a handful of cases. In Fig. [3] 
we indicate the results found for tw o bright, lensed LBGs at 
z ~ 3 from lSiana etaT] ([2008, 2009). These objects have di- 
rect detections at (observed) 24 and 70/im (as well as 850/im 
and 1.2 mm in the case of cB58), allowing a good estimate 
of Ljr. It has been noted that both objects appear to lie sub- 
stantially below IRXm99,o> perhaps requiring a modified ex- 
tinction curve or a n extremely you ng stellar population (see 
iBaker et ail 120011 : ISiana et ail 120091 and the discussion on 
"young galaxies" by R10). An extinction curve more simi- 
lar to that of the Small Magellanic Cloud (SMC), i.e., extinc- 
tion that occurs in a more uniform rather than a patchy fore- 
ground dust screen would require less emission by dust for 
the same stellar population age, and hence lower IRX (see red 
curve in Fig. [J). cB58 lies on the lower envelope formed by 
LBAs and z ~ 2 galaxies, and the observed scatter is larger 
than the difference expected between the standard Calzetti- 
type and SMC-type extinction curves for small j3. The cosmic 
Eye, however, has a much larger slope (/3 » —0.3) and lies 
well below the starburst IRX relations. Previous studies have 
shown that offsets from IRX are comm only due to mor e com- 
plicated star formation histories (e.g. Kong et al. 2004). This 
is a likely explanation, at least for t he local SFGs (left p anel), 
BzK galaxies (right panel, see also iNordon et al.ll2010l) . and 



LBAs (see objects marked "DCOs" in Fig. [Hand O09) hav- 
ing relatively low IRX ratios. For the cosmic Eye, however, a 
modified extinction curve may indeed be required in order to 
simultaneously explai n its relatively sim ple stellar population 
and its low IRX ratio (ISiana et alJl2009l) . 

4. DISCUSSION AND CONCLUSIONS 

The first goal of this paper was to investigate whether LBAs 
are consistent with the M99 relation, and we have shown that 
this is indeed the case. This is an important result: Because 
LBAs are much better analogs to LBGs compared to previous 
local samples, including M99, the applicability of the local 
dust correction at high redshift now appears to be on much 
more solid ground. Confirmation of this result is given by the 
fact that LBAs are also very similar to typical LBGs in terms 
of their dust properties (Figs. |2]and[3]), although some caution 
must be taken as true Lm measurements are still sparse at 
high redshift. 

We can show the implications of these results on the deter- 
mination of the SFR density at z ~ 3 — 6. B09 used large 
dropout samples from the Hubble Ultra Deep Field to quan- 
tify the relation between pw and Mjjv, and computed the to- 
tal dust corrections to the SFR density by integrating the UV 
luminosity functions a t each redshift (Bouwens et al J 120071 : 
iReddy & S teidel 2009) and assuming that IRXm99.o is valid. 
The method and data that we use here is similar to that of B09 
(Fig. 8 and §5.4 of that work), but instead of using IRXm99.o 
we show the results for our new IRX relations derived in §3. 1 . 
The results are shown in Fig. |4j which gives the total correc- 
tion factor that is needed to dust-correct a measurement of the 
UV luminosity (or SFR) density at z ~ 3 — 6. Not surpris- 
ingly, IKXlba (thick black lines) and IRX^/gg, inner (dotted 
lines) give quite similar results to the relation used by B09 
(thin dashed lines). Note that if we assume that IRXm99, total 
measured within the larger aperture provides a better estimate 
of the total attenuation, the total dust-correction would have to 
be lowered by ~ 0.2 dex at all redshifts and for all luminosi- 
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Fig. 3. — : IRX-/3 at low (left) and high redshift (right). LBAs are indicated in both panels ( black points). Left panel: M99 
starbur sts with R < R mnx (grey squares), SFGs (cyan upside-down triangles; Bua t et al.l ( 120101) ). and (U)LIRGs (open magenta 
circles: lHowell et aTl (120101) . Right panel: LBGs (filled blue triangles), BzK galaxies (open purple triangles), distant red galaxies 
(open red squa res) and SMGs (open or ange pentagons) at z ~ 2 — 3 (Reddy et al. 2006. R10), lensed LBGs at z ~ 3 (large 
green triangles: ISiana et al.l (|2008, 2009), and statistical detections of LBGs at z ~ 3 detected at 1 .1mm (lower magenta triangle, 
Magd iTet al.l l2010ah or at 160 /urn with Herschel (upper magenta triangle, Magd is"et al.ll20 10b). Curves show IRXm99. inner 
(dotted), IRXa/99, total (dashed) and IKXlba (thick black) from Fig. [Hand §3.1. The IRX relation assuming a SMC extinction 
curve is indicated by the thin red line. 
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Fig. 4. — : The multiplication factors needed to dust-correct the SFR density at z ~ 3 — 6 based on the UV slopes and luminosity 
functions (LF) of U, B, V, i "dropout" galaxies from B09, and assuming that the local attenuation corrections are valid at high 
redshift. Panels shows the total dust corrections implied when integrating the LFs over (L/L* =3 )>0.3, 0.04, and 0. The curves 
show the results when using IRX^s^ (black solid line), IRXM99,inner (g re y dotted line) and IRXm99, totai (grey dashed line) 
from §3.1. The thin dashed curve shows the results of B09 using the original M99 relation (IRXj\/gg o). 



ties (grey dashed lines). We conclude that current determina- 
tions of the cosmic SFR density based on dust-corrected UV 
data from Lyman break galaxy surveys appear to be on solid 
ground, but, as noted by B09, this is mainly due to the large 
contribution from very faint, blue sources at high redshift. 

Direct measurements of Lm for relatively red {(3 > 0) 
galaxies at z > 2 — 3 are needed to constrain the IRX rela- 
tion in this regime, while the most obscured sources typically 
do not follow the IRX relations at all (e.g. ULIRGs at low 



and SMGs at high redshift). This illustrates the importance 
of surveys in the IR currently being performed by Herschel. 
In forthcoming papers we will present a more detailed anal- 
ysis of the dust properties of LBAs and LBGs based on IRS 
spectra as well as new far-IR data from Herschel. 

We thank Guinevere Kauffmann, Rychard Bouwens, Luca 
Cortese and Ranga-Ram Chary for providing useful com- 
ments. 
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TABLE 1 

Lyman Break Analogs: UV and IR Measurements. 
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ID rdcrs lo the labels shown in Figure 1. 
' Far-IR luminosity integrated over the range 40- 120 ^im. 

IR luminosity integrated over the range 3-1000 ^im. 
' Indicates objects referred to as Dominant Compact Objects in O09. 
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